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The computational complexity of time-dependent perturbation theory is well-known to be largely 
combinatorial whatever the chosen expansion method and family of parameters (combinatorial se- 
£SJ ■ quences, Goldstone and other Feynman-type diagrams...). We show that a very efficient perturba- 

tive expansion, both for theoretical and numerical purposes, can be obtained through an original 
parametrization by trees and generalized iterated integrals. We emphasize above all the simplicity 
and naturality of the new approach that links perturbation theory with classical and recent results 
in enumerative and algebraic combinatorics. These tools are applied to the adiabatic approximation 
and the effective Hamiltonian. We prove perturbatively and non-perturbatively the convergence of 
Morita's generalization of the Gell-Mann and Low wavefunction. We show that summing all the 
terms associated to the same tree leads to an utter simplification where the sum is simpler than any 
of its terms. Finally, we recover the time-independent equation for the wave operator and we give 
an explicit non-recursive expression for the term corresponding to an arbitrary tree. 
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. T. INTRODUCTION 
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Effective Hamiltonians provide a way to determine the low-energy eigenvalues of a (possibly infinite dimensional) 
Hamiltonian by diagonalizing a matrix defined in a subspace of small dimension, called the model space and hereafter 
denoted by M. Because of this appealing feature, effective Hamiltonians are used in nuclear, atomic, molecular, 
5^ ' chemical and solid-state physical. 

These theories are plagued with a tremendous combinatorial complexity because of the presence of folded diagrams 
(to avoid singularities of the adiabatic limit), partial resummations, subtle "linkedness" properties and the exponential 
\ growth of the number of graphs with the order of perturbation. This complexity has two consequences: on the one 
hand, few results are proved in the mathematical sense of the word, on the other hand, it is difficult to see what is 
the underlying structure of the perturbative expansion that could lead to useful resummations and non-perturbative 
approximations. 

To avoid these pitfalls, we take a bird's-eye view of the problem and consider a general time-dependent Hamiltonian 
H(t). This way, we disentangle the problem from the various particular forms that can be given to the Hamiltonian 
and which have lead in the past to various perturbative expansions. To be precise, take the example of fcrmions 
in molecular systems. The Coulomb interaction between the electrons (say V) can be viewed as a perturbation 
of a "free Hamiltonian" modeling the interaction with the nuclei (in the Born-Oppenheimer approximation). One 
can take advantage of the particular form of V (which is a linear combination of products of two creation and two 
annihilation operators in the second quantization picture) to represent the perturbative expansions using a given family 
of Goldstone diagrams (see e.g. ref. Hj for such a family). However, the general results on perturbative expansions 
(such as the convergence of the time-dependent wave operator) do not depend on such a particular choice. 

Thus, we consider an Hamiltonian H(t) and we build its evolution operator U(t,to), which is the solution of the 
Schrodinger equation (in units h = 1) 

= H(t)U(t,t ), (1) 

with the boundary condition Ufa, to) = 1. In perturbation theory, H(t) := e _e l*le zi?ot ye _lffot is the adiabatically 
switched interaction Hamiltonian in the interaction picture (here Hq and V stand respectively for the "free" and 
interaction terms of the initial Hamiltonian) and singularities show up in the adiabatic limit (to — > — 00 and e — > 0). 
Morita discovered? that, in this setting, the time- dependent wave operator 

n(t,t ) := u(t,t )p(pu(t,t a )py\ 



in 



>< 



where P is the projection onto the model space M, has no singularity in the adiabatic limit. Moreover, the wave 
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operator determines the effective Hamiltonian because 



H eS := limPFfi(0,-oo). 



e->0 



(2) 



However, as we have aiready alluded to, the effective computation of these operators raises several combinatorial and 
analytical problems that have been addressed in a long series of articles (several of which will be referred to in the 
present article). 

In the first sections of the paper, we consider a general time-dependent Hamiltonian H(t) (not necessarily in the 
interaction picture) . In this broader setting, Jolicard 3 - found that the time-dependent wave operator provides also a 
powerful description of the evolution of quantum systems (see ref.0 for applications). Then, we derive three (rigorously 
proven) series expansions of the wave operator. The first one is classical and can be physically interpreted as the 
replacement of causality (i.e. the Heaviside step function 6(t — t')) by a "propagator" 6 P (t — t') := 9(t — t')Q — 8(t' — t)P, 
where Q = 1 — P. This "propagator" is causal out of the model space (0p(t — t') = for i < i' on the image of Q) and 
anticausal on it, like the Feynman propagator of quantum field theory^. However, this sum of causal and anticausal 
orderings is cumbersome to use in practice. A second series expansion is obtained by writing the wave operator 
as a sum of integrals over all possible time orderings of the Hamiltonians H(ti) (see Sect. IIIB[) . This expansion, 
parametrized by all the permutations (or equivalent families), is used in many-body theory and gives rise to a large 
number of complicated terms. The third expansion is obtained by noticing that some time orderings can be added to 
give simpler expressions. This series is naturally indexed by trees and is the main new tool developed in the present 
paper. Among others, we derive a very simple recurrence relation for the terms of the series. We also show that the 
very structure of the corresponding generalized iterated integrals showing up in the expansion is interesting on its 
own. These integrals carry naturally a rich algebraic structure that is connected to several recent results in the field 
of combinatorial Hopf algebras and noncommutative symmetric functions. The corresponding algebraic results that 
point out in the direction of the existence of a specific Lie theory for effective Hamiltonians (generalizing the usual 
Lie theory) are gathered in an Appendix. 

In the last sections of the paper, we restrict H(t) to the interaction picture and we consider the adiabatic limit. 
We first prove that the adiabatic limit exists non perturbatively. We show that the effective Hamiltonian defined 
by eq. ([5]) has the expected properties. Then, we expand the series and we give a rigorous (but lengthy) proof that 
the term corresponding to each time ordering has an adiabatic limit. Then, we consider the series indexed by trees 
and we give a short and easy proof of the existence of that limit. Finally, we provide a direct rule to calculate 
the term corresponding to a given tree and establish the connection between the time-dependent approach and the 
time-independent equations discovered by LindgrenS. and Kvasnickai. 

The existence of this series indexed by trees can be useful in many ways: (i) It describes a sort of superstructure 
that is common to all many-body theories without knowing the exact form of the interaction Hamiltonian; (ii) 
It considerably simplifies the manipulation of the general term of the series by providing a powerful recurrence 
relation; (iii) It provides simple algorithms to calculate the terms of the series; (iv) The number of trees of order n, 



where the Hamiltonian Hq is not quadratic. Indeed, many-body theories most often require the Hamiltonian Hq to 
be free, i.e. to be a quadratic function of the fields 1 -. As noticed by Bulaevskii 1 — , this is not the good point of view 
for some applications. For example, in the microscopic theory of the interaction of radiation with matter, it is natural 
to take for Ho the Hamiltonian describing electrons and nuclei in Coulomb interaction^, the perturbation being the 
interaction with the transverse electric field. In that case, quadratic free Hamiltonians many-body theories break 
down whereas our approach is still valid. Actually, it is precisely for that reason that we originally developed the tree 
series approach; (vi) Last, but not least, the tree-theoretical approach connects many-body theories with a large field 
of knowledge that originates in the "birth of modern bijective combinatorics" in the seventies with in particular the 
seminal works of Foata, Schiitzenberger and Vienno t 13 ' 14 . See e.g. ref. Qj5l for a survey of the modern combinatorial 
theory of tree-like structures. 

From the physical point of view, the tree expansion is particularly interesting in the adiabatic limit. Indeed, the 
denominator of each of its terms is a product of Ef — factors, where E^ is the energy of a state in the model 
space M and the energy of a state not belonging to M. In the usual many-body expansions, the denominators are 
products of E® — ^ . E^ factors, where the sums contain various numbers of elements (corresponding therefore 
to multiple transitions between low-energy and cxcitated levels). In that sense, the tree expansion is the simplest 
possible because each term is a product of single transitions between two states. 

We now list the main new results of this paper: (i) A recursion formula that generates the simplified terms of 
the time-dependent perturbation series (theorem ; (ii) when the interaction is adiabatically switched on, a non- 
perturbative proof of the convergence of the wave operator and a characterization of the states of the model space 
that are transformed into eigenstates of H by the wave operator (theorem O; (iii) a proof of the existence of the 




being subexponential, it improves the convergence of the series^; (v) It can deal with problems 
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adiabatic limit for the terms of the series expansion of the wave operator (theorem [9]); (iv) a recursive formula 
f lemma [T3"f and an explicit form (theorem I14[) for the general term of the time- independent perturbation series. 

II. TIME-DEPENDENT HAMILTONIAN AND COMBINATORICS 

We consider a time-dependent Hamiltonian H(t), which is a self-adjoint operator on a Hilbert space Ti, and its 
evolution operator U(t,to) defined in eq. ([1}. Since we are interested in the combinatorial aspects of the problem, we 
consider the simple case where H (t) is a strongly continuous map from R into the bounded self-adjoint operators on 
H— . In that case, the Picard-Dyson series 

00 pt /^l ptn — 1 

1 + / d< i / d< 2---/ dt n H(t 1 )...H(t n ), 

n — i J to "to ^t 

converges in the uniform operator topology to U (t, to) and U(t,to) is a jointly strongly continuous two-parameter 
family of unitarics on % (see section X.12 of ref. \Ym . 

Following Morita£, Jolicard^ established a connection between the evolution operator and the effective Hamiltonian 
approach by defining 

n(t,t ) := u(t,t )p(pu(t,t )py 1 , 

where P is a projection operator onto M and (PU(t,to)P) 1 is the inverse of PU(t,t )P as a map from M = PH 
to itself. This map is invertible if and only if there is no state \<j>) in M such that (<j>\U(t,to)P — 0. This condition is 
similar to the one of time-independent perturbation theory^. We assume from now on that the condition is satisfied 
and we define three expansions for 0(t, to)- 

A. First expression for Q 

We start by proving an elegant expression for J7, that was stated by Michels and Suttorpi^ and Dmitriev and 
Solnyshkina^. 

Theorem 1 

n(t,t ) = P + Qj2(-t) n dh dh... dt n H(t 1 )9 P (t 1 ~t 2 )H(t 2 )...8 P (t n _ 1 -t n )H(t n )P, (3) 
n= i Jto "to "to 

where Q = 1 — P and 0p(t) = 9(t)Q — 6{—i)P, with 9 the Heaviside step function. 

Proof. We first rewrite the Picard-Dyson series as U(t,to) = 1 + Yl n Un(t,to) with U\(t,to) :— —i J t dtiH(t\) and, 
for n > 1, 

U n (t,t ) := (-*)" f dh... f dt n H(h)6(h-t 2 )...8(t n _ 1 -t n )H(t n ). 
Jto Jto 

Then, by using 8(t) + 0(—i) = 1, wc notice that 

9{t) = P + 9(t) -9(t)P -9(-t)P = P + 9 P (t). 

Now, we replace 9(t) by the sum of operators P + 6p(t) in the expression for U n (t, t ). This gives us 2 n_1 terms with 
various numbers of P and 9p. Denote by C n (t, to) the term with no P (with the particular case Ci(t, to) = Ui(t, to)). 
Take then any other term. There is an index i such that the first P from the left occurs after H(h). Therefore, the 
integrand of this term is 

H(h)9 P (h -h)... 9 P (t i _ 1 - h)H(t l )PH(t l+1 ) . . . 

Observe that the integral over t%, . . . ,t% is independent from the integral over tj+i, . . . ,t n . The first integral gives 
Ci{t,to), the second integral is a term of the Picard-Dyson series for U n _i(t,t ). Thus, the sum of the 2™ _1 terms 
yields 

n-l 

U n (t,t ) = C„(i,io) + ^2a(t,to)PU n - l {t,to). 
i=l 
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If we denote by K(t, to) the sum of all the C n (t, to) with n > 0, we obtain U = 1 + K + KP(U — 1), so that 

UP = P + KPUP. (4) 

The operator K is called the reduced evolution operator by Lindgren and collaborators 2 —. If we define u> := P + QKP, 
eq. (@| becomes 

UP = P + {uj- P + PKP)PUP = P + ujPUP - PUP + PKPUP. 
This equation can be simplified by using eq. (|4]) again 

UP = P + ujPUP — PP = ujPU P. 
Thus, lu = fi and eq. ([3]) is satisfied. □ 
Despite its elegance, eq. (j3]) is not immediately usable. To illustrate this point, consider the third-order term 

n 3 = iQ f dt x f dt 2 f dt 3 H{ti)0p(t 1 -t2)H(t 2 )9p(t2-t 3 )H(t 3 )P, 

Jto Jto Jto 

If we expand 9p(t) = 6(t)Q — 6(—t)P, we obtain four terms 

tQH(ti)QH(t 2 )QH(t 3 )P for t x > t 2 and t 2 > t 3 , 

-iQH{ti)QH{t 2 )PH(t 3 )P for t x > t 2 and t 2 < t 3 , 

-iQH{ti)PH(t 2 )QH(t 3 )P for h < t 2 and t 2 > t 3 , 

iQH(tx)PH{t 2 )PH{t 3 )P for t x < t 2 and t 2 < t 3 . 

The first and last terms have integration range t\ > t 2 > t 3 and ^3 > t 2 > tx, respectively and give rise to iterated 
integrals. The integration range of the second term is t\ > t 2 and t 2 < t 3 . Such an integration range is not convenient 
because the relative position of t\ and £3 is not specified. The integration range has to be split into the two subranges 
ti > £3 > t 2 and t 3 > t\ > t 2 . Each subrange defines now an iterated integral. For example t\ > t 3 > t 2 gives 

-1 f dh f 1 dt 3 f 3 dt 2 QH{h)QH{t 2 )PH(t 3 )P. 

Jto Jto Jto 

Similarly, the integration range of the third term {t\ < t 2 and t 2 > t 3 ) is the union of t 2 > t\ > t 3 and t 2 > t 3 > t\. 
We see that ^3 is sum of six iterated integrals corresponding to the six possible orderings of ti, t 2 and t 3 . 



B. Q in terms of permutations 

We consider again the previous example, and we change variables to have a fixed integration range s\ > s 2 > s 3 . 
If we sum over all time orderings, we obtain 

CL 3 = 1 f dsi f 1 ds 2 I 2 ds 3 (QH(s 1 )QH(s 2 )QH(s 3 )P-QH(s 1 )QH(s 3 )PH(s 2 )P-QH(s 2 )QH(s 3 )PH(s 1 )P 

Jto Jto Jto 

-QH(s 2 )PH( Sl )QH(s 3 )P - QH{s 3 )PH{s x )QH{s 2 )P + QH{s 3 )PH{s 2 )PH{si)P) . (5) 

In ref. we showed that this result can be generalized to all orders and that £l n is a sum of nl iterated integrals 
corresponding to all the orderings of t%, . . . , t n . More precisely, we obtained the series expansion for 51 

fi(i,*o) = P + J2 E / f 1 dt 2 ... I*" 1 dt n X{t a{1) )...X{t a(n) ), (6) 

where S n is the group of permutations of n elements. The operators X are defined, for n = 1, by X(t) := —iQH{t)P 
and, for n > 1 and any a e S n , by X{t,y^) :— —iQH^u^) and 

X(t a ( p )) := -iQH(t a{p) ) if 1 < p < n and a(p) > a(p - 1), 

X(t a ( p )) := tPH(t a ( p )) if I < p < n and a(p) < a(p - 1), 

X(t a(n) ) := -iQH(t a(n) )P ii a{n) ><r(n - 1), 

X{t a{n) ) := iPH(t a(n) )P if <r(n) < <r(n - 1). 
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Each term of eq. (0 is now written as an iterated integral. However, the expansion © is still not optimal because 
some of its terms can be summed to get a simpler expression. As an example, consider the fourth and fifth terms of 
eq. ([5]), where we replace s, by h. 

Z := -i f dh f 1 dt 2 f * dt 3 (QH(h)PH(h)QH(h)P + QH(h)PH(h)QH(h)P)- 

J to J to J to 

The first and second terms of the right hand side of this equation are denoted by Z\ and Z 2 , respectively. We 
transform Z 2 by exchanging variables h and t 3 . 

Z 2 = -1 [ dh [ dt 3 I* dt 2 QH{h)PH{h)QH{h)P = -1 f dh [ * dt 2 j * dt 3 QH(t 2 )PH(h)QH(t 3 )P, 

J to <J to J to J to <J to J £2 

where we also exchanged the order of the integrations over t 3 and t 2 . This can be added to Z\ and we obtain the 
simpler expression 

Z = -i f dh f 1 dt 2 f 1 dt 3 QH(t 2 )PH(h)QH{h)P. 

J to Jto J to 

Such a simplification is not possible for the other terms of Q 3 . In the next section, we determine how this simplification 
can be extended to the general term of 57. 

Before closing this section, we need to specify more precisely the relation between the permutations a and the 
sequence of P and Q in the expansion of eq. ([3]). When we expand all the Qp(ti — i^+i) in eq. ([3]), we obtain an 
integrand of the form (—i) n QH{h)RiH(t 2 ) . . . R n -iH(t n )P, multiplied by a product of Heaviside functions, where 
Ri takes the value —P or Q. We aim to determine the relation between the sequence R\ . . . R n -i and the permutations 
a in eq. ©. From the definition of X{t a (^y), it appears that Ri = —P if a{i) > a(i + I) and Ri — Q if a{i) < a(i + 1). 
The set of indices i such that a(i) > a(i + 1) is called the descent set of <r, denoted by D a . It is also called the 
shape of the permutation^ 3 .. For instance, the descent set of the permutations (213) and (312) is {1}, corresponding 
to(Ri,R 2 ) = (-P,Q). 

C. Permutations and trees 

In many-body physics, the expansion in Goldstone diagrams corresponds (among other things) to the expansion of 
fl into all time orderings In that context, several authors noticed that some diagrams corresponding to different 
orderings can be added to give a simple sum^ ' 19 ' 24 — , as we saw in the previous section. These are special cases of the 
simplification that we shall present which, as far as we know, was never stated in full generality. The first difficulty 
is to find the proper combinatorial object to represent the sets of permutations that lead to simplified sums. We 
shall find it after a short algebraic detour meant to recall the notion of tree and its relation to permutation s 13 ' 23 . 
The trees we consider have received various names in the literature: binary trees in ref. |23| . but also (to quote only a 
few occurrences) plane rooted complete binary tree 2 ^, extended binary trees 2 ^ or planar binary freest. Since these 
objects are rarely used in physics, we first link them with the trees commonly found in graph theory. This lengthy 
definition will then be replaced by a much easier one. 

A common tree is a connected graph without loops. In other words, a common tree is a set of vertices linked by 
edges such that there is a unique way to go from any vertex to any another one by travelling along edges. An example 
of a common tree is given in fig. [TJa). A rooted tree is a common tree where one vertex was selected. This particular 
vertex is called the root. The level of a vertex in a rooted tree is the number of edges that separates this vertex from 
the root. The root has level 0. It is natural to draw a rooted tree by putting the root at the bottom and drawing a 
vertex above another one if its level is larger. The rooted tree of fig. [D(b) was obtained by selecting as the root the 
lowest left vertex of fig. [TJa). The root is indicated by a white dot and a dangling line. In a rooted tree, the children 
of a vertex v are the vertices v' that are linked to v by an edge and such that the level of v' is larger than the level 
of v. A plane rooted binary tree is a rooted tree where each vertex has zero, one or two children, and each edge is 
oriented either to the left or to the right. If a vertex has two children, then one of the edges is oriented to the left 
and the other one to the right. Fig. [IJc) shows one of the plane rooted trees that can be obtained from the rooted 
tree of fig- [2(b). The adjective "plane" means that an edge going to the right cannot be identified with, or deformed 
into an edge going to the left. A plane rooted complete binary tree is a plane rooted binary tree where each vertex is 
"completed" by drawing leaves as follows: if a vertex has no child, draw a leaf (i.e. a line) to the left and one to the 
right, if a vertex has one child, then draw a leaf to the right if the child is to the left and draw it to the left if the 
child is to the right. Fig. HJd) shows the plane rooted complete binary tree that is obtained from the plane rooted 
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(a) (b) (c) (d) (f) 

FIG. 1: Construction of plane rooted complete binary tree, (a) a common tree; (b) a rooted tree with its vertex levels; (c) a 
plane rooted binary tree; (d) a plane rooted complete binary tree; (f) a simplified plane rooted complete binary tree. 

binary tree of fig. HJc). In practice, the vertices are no longer necessary and they are not drawn, as in fig. [TJf ) . If 
Y n denotes the set of plane rooted complete binary trees with n vertices, we see that Y\ = { Y }, Y 2 = { I , I }, 

Y3 = { J , [ , ^Y' , T*' i Y^ }• They are much more numerous than the common trees (there is only one common 
tree with one, two or three vertices). For notational convenience, plane rooted complete binary trees will be simply 
called "trees" in the rest of the paper. 

Fortunately, there exists a much simpler definition of the trees, that use a sort of building rule. We first denote 
the empty tree, i.e. the tree with no vertex, by | , which is a dangling line without root (a dangling line with a root 

and no other vertex belongs to the tree Y )• Then, for any integer n, Y n , is defined recursively by Yq := { I } and, for 
n > 0, Y n := {Ti V T 2 : T x E Y kl T 2 E Y n _ k _i, k = 0, . . . , n - 1}, where T x V T 2 is the grafting of the two trees T x 
and T 2 , by which the dangling lines of T± and T 2 are brought together and a new root (with its own dangling line) is 

grown at their juncture. For example, | V ] = Y,|vY = f,TvY = ""V ■ Note that each tree of Y n has n 
vertices (including a root) and n + 1 leaves. The order \T\ of a tree T is the number of its vertices. 

If C n denotes the number of elements of Y n , the recursive definition of Y n implies that Co = 1 and C n — 
Sfe=o CkCn-k-i, so that C n — ^y( 2 ™) are the famous Catalan numbers. For n=0 to 10, C„ = l, 1, 2, 5, 14, 
42, 132, 429, 1430, 4862, 16796. The Catalan numbers enumerate a large number of (sometimes quite different) 
combinatorial objects 3 -!. The main practical interest of trees with respect to other combinatorial interpretations is 
that their recursive definition is very easy to implement. 

We noticed in the previous section that, at order three, the terms fl a corresponding to two specific permutations a 
can be added to give a simple result. At the general order, we shall see that the sum of f2 CT is simple if it is carried 
out over permutations a associated with the same tree. But for this, we need to associate a tree to each permutation. 
The relevant map from permutations to trees belongs to "the ABC's of classical enumeration"— and is historically 
one of the founding tools of modern bijective and algebraic combinatoric o 13 ' 23 ' 33 . We describe it in the next section. 



1. From permutations to trees 

We first map any n-tuple I = . . . , i n ) of distinct integers to a tree T E Y n . The mapping (j) is defined recursively 
as follows: if I = (i) contains a single integer, then (f)(1) := Y ; otherwise, pick up the smallest element ik of /, then 

0(7) := I V0((i 2> ...,i n )), if k = l, 

(f>(I) := 0((ii,...,i n _i)) V I , if k = n, 

(f>(I) := . . ,z fc _i)) V <f)((i k+ i,...,i n )), otherwise. 
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In the following, we frequently abuse notation by writing 0(«i, . . . , i n ) or even 0(ii ...i n ) instead of 0((«i, ■ ■ . , i n )) ■ 
For a permutation a G S n , the corresponding tree is 0(c) = 0(cr(l), . . . , cr(n)J. For example 



0(1) = 


Y, 






0(12) = 


r, 


0(21) = ^ , 




0(123) = 




0(132)= f , 


0(213) = 


0(231) = 




0(312) = Y . 


0(321) = 



Note that the two permutations (213) and (312) correspond to the same tree, and they are also the two permutations 
that add up to a simple sum in the calculation of Z at the end of section Hi Bl This is not a coincidence. 

To simplify the proofs, we embrace the three cases of the definition of into a single one as follows. We first define 
the concatenation product of two tuples I = (ii, . . . , i n ) and J = (ji, . . . ,j m ) as / • J = (ii, . . . ,i n , ji, . . . , j m ). We 
extend this definition to the case of the zero-tuple Io = by Iq ■ I = I ■ Io = I ■ Then, for any n-tuple / of distinct 
integers, we define (f)(1) by 0(7) — | if n = and 0(7) = (f)(I\) V (f)(Ii) otherwise, where I\ and Ii are determined by 
I = I\ ■ (mini) • Ii. Note that I\ or Ii can be the zero-tuple. 

We first prove an easy lemma. 

Lemma 2 If the elements of the two n-tuples I = (ii, . . . , i n ) and J — (ji, . • . ,j n ) of distinct integers have the same 
ordering (i.e. if ik < ii if and only if jk < ji for all k and I in {1, . . . ,n}), then 0(1) = 4>(J). 

Proof. The proof is by induction. If n = 0, then I = J = and 0(7) = 0(J) = I • Assume that the property 
is true for fc-tuples of distinct integers up to k = n — 1 and take two n-tuples / and J having the same ordering. 
Then, the minimum element of both is at the same position k (i.e. mini = ik and min J = jk) and I = I\ ■ (ik) ■ Ii, 
J = J\ ■ (jk) ■ Ji, where I\ and J\ (Ii and J2, respectively) are two (k — l)-tuples ((n — fc)-tuples, respectively) of 
distinct integers have the same ordering. By the recursion hypothesis, we have 0(/i) = 0(Ji) and <f>(Ii) — <t>(Ji) and 
the definition of gives us (f)(1) = 0(ii) V 0(I 2 ) = 0(J X ) V 0(J 2 ) = 0(J). □ 

As a useful particular case, we consider the situation where J describes the ordering of the elements of I: we 
order the elements of / = (i±, . . . ,«„) increasingly as ik t < ■ ■ ■ < ik„- Then ji is the position of i\ in this ordering. 
More formally, J := (r _1 (l), . . . , r _1 (n)), where r is the permutation (k±, . . . , k n ). The n-tuple J is called the 
standardization of / and it is denoted by st(I). If we take the example of / = (5, 8, 2), the position of 5, 8 and 2 in 
the ordering 2 < 5 < 8 is 2, 3 and 1, respectively. Thus, st(5, 8, 2) = (2, 3, 1). By construction, / and st(I) have the 
same ordering and (f)(1) — (f>(st(I)). We extend the standardization to the case of / = by si(0) = 0. 



2. From trees to permutations 

Conversely, we shall need to know the permutations corresponding to a given tree: St '■= {c £ S\t\ ■ 0( cr ) = T} 

(we extend this definition to the case of T = \ by defining the zero-element permutation group Sq := {0}). The 
solution of this problem is given by 

Lemma 3 If T = T\\J Ti, where \T\\ = n and \Ti\ — m (n or m can be zero), all the permutations of St have the 
form I = I\ ■ (1) • Ii, where I\ is a subset of n elements of {2, . . . , n + m + 1} ordered according to a permutation a 
of St x (i.e. st(Ii) — a) and Ii is the complement of I\ in {2, . . . , n + m + 1}, ordered according to a permutation /3 
of St 2 (i.e. st(Ii) — (3). 

Proof. The proof is given in refs. and[3(J but we can sketch it here for completeness. The simplest examples are 
S T = {0} for T = \ and S T = {(1)} for T = Y ■ Now take T = T x V T 2 as in the lemma. By the definition of 
0, the minimum of the tuple / = (c(l), . . . ,o(n + m + 1)) is a(n + 1) = 1 and I = 1% • (1) • Ii, where 0(/i) = T\ 
and 4>(Ii) — Ti. We saw in the previous section that 4>(I\) — 4>(st(I\)). By definition st(I\) is a permutation of S n . 
Therefore, st(I\) belongs to Stx and, similarly, st(Ii) belongs to St 2 - It is now enough to check that each element of 
St is obtained exactly once by running the construction over all orderings and all permutations of T\ and Ti. □ 

This lemma allows us to recursively determine the number of elements of St, denoted by |«St|, by \St\ = 1 for T = 
and T — Y and, for T = T x V T 2 , 

\St\ = O^'^IStAIStJ. (7) 
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See ref. Hjfij for an alternative approach. 

Example: Consider the tree T = Ti V T 2 , with T\ = Y and T 2 = ■ so that n — 1 and m = 3. St x contains 
the single permutation a = (1) and, according to the examples given in the previous section, the two permutations 
of St 2 are fli = (213) and j3 2 = (312). We choose the permutations a and we pick up n — 1 element (for example 
3) in the set J = {2, 3, 4, 5}, so that I\ = (3) and we order the remaining elements {2, 4, 5} according to /3i, so that 
I 2 = (4, 2, 5). This gives us a — (31425). If we pick up the other elements of J to build I\ we obtain (21435), (41325) 
and (51324). We add the elements obtained by choosing f3 2 and we obtain eight permutations: 

S T = {(21435), (21534), (31425), (31524), (41325), (41523), (51324), (51423)}. 

We can check that eq. Q holds and that \St\ = 8. 

D. Recursion formula 

The permutations corresponding to a tree can be used to make a partial summation of the terms of the Picard-Dyson 
expansion. 

Definition 4 For any tree T, we define D,T(t,to) by r2r(i,to) —PifT— | and 

dti / d; 2 .../ dt n x{t a{l) )...x{t a{n) ), 

otherwise, where n — \T\. 

With this notation we have obviously ft(t,to) — ^ T f2^(t, t ) = J2a^^^^o)j with the notation fl a (t,t ) := 
J t * dti dt 2 ■ ■ ■ J*™ 1 dt n X{t a ^) . . . X(t a ( n )) and where a runs over all permutations (of all orders). The term 
of order of this series is fi| = P and the term of order one is 



nr(*,to) = - 




dsQH(s)P, 



for T = Y ■ The other terms enjoy a remarkably simple recurrence relation: 
Theorem 5 // \T\ > 1, then fiT(i,io) can be expressed recursively by 

n T (t,t ) = -1 [ dsQH{s)n T2 {s,t ), if T = I V T 2 , 
Cl T (t,t ) = 1 f dsn Tl (s,t Q )H(s)P, if T = T 1 y 1 , 

Jto 

n T (t,t ) = % [ d s n Tl (s,t )H(s)n T2 ( s ,t Q ), xt = t 1 \/t 2 , 

where Ti ^ \ and T 2 ^ \ . 

Note that a similar recursive expression was conjectured by Olszewski for the nondegenerate Rayleigh-Schrodinger 
expansion^!. 

Proof. Let us prove the theorem recursively. Consider an arbitrary T = T\ V T 2l |T| > 1, and assume the formulas to 

hold for all the trees T' with |T'| < |T|. Consider for example the case where T\ ^ \ and T 2 ^ \ (the other cases 
are even simpler). We define: 

A T :=ij dsn Tl {s,t a )H(s)n T2 (s,t ) = ^ 1 f dsCt a (s,t )H(s)np(s,t ). 
•' t ° aes Tl ,0eS T2 

This first important point is that, for a given tree T, all the permutations a 6 St have the same descent set. This is a 
well-known fac t 23 i 30 that can be deduced from the characterization of St at the end of section lll CI As a consequence, 
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the sequence of operators P, Q and H is the same for all a and (3 in At, and only the order of the arguments t, 
varies. Therefore, the conditions of lemma [T51 fsee appendix IA"| are satisfied and we get: 

4r = ^fi 7 (t,t ), 

7 

where 7 runs over the permutations such that 7(|T\ + 1) = 1, st(j(l), ...,7(|Ti|)) £ S^, st(j(\Ti\ + 2), ...,7(|Ti| + 
\T 2 \ + 1)) £ St 2 - The set of permutations 7 satisfying these equations is precisely St, so that, finally: A T = ^t(^ ^o)- 
This concludes the proof of the theorem. □ 



E. Remarks 

1. Nonlinear integral equation 
If we denote x(^o) = ^(ijio) — P> then the recurrence relations add up to 

» x (t,to) - / dsQH(s)P+ f dsQH(s) X (s,t ) - ( ds X (s,t )H(s)P - ( ds X (s,t )H(s)x(s,t ). (8) 

J to "^0 J to J to 

The derivative of this equation with respect to t was obtained in a different way by Jolicard 3 . 

2. Permutations, trees and descents 

We saw that, for a given tree T, all permutations of St have the same descent set. We can now give more detail o 30 ' 34 . 
The relation between the trees and the sequences of operators P and Q in eq. ([3]) is very simple. Consider the sequence 
of leaves from left to right. Each leaf pointing to the right corresponds to a P, each leaf pointing to the left correspond 

to a Q. For example, the tree ] corresponds to the sequence QQPP. From the combinatorial point of view, this 
description emphasizes the existence of a relationship between trees and descent sets (or, equivalently, hypercubes), 
see e.g. refs. [23ll34| - [37l and our Appendix. 

III. ADIABATIC SWITCHING 

Morita's formula is most often applied to an interaction Hamiltonian H e (t) :— Q~' i \ t \Q lH o t Ve~ lHot . We write 
Eq, E n , ... and $0, ■•■ for the eigenvalues of Hq and for an orthogonal basis of corresponding eigenstates. We 

assume that the spectrum is discrete and that the eigenvalues are ordered by (weakly) increasing order. The ground 
state may be degenerate (E n = E\ = ... = Ef. for a given k). The model space M (see the Introduction) is the vector 
space generated by the lowest TV eigenstates of H (with N > k). We assume that the energies of the eigenstates of 
M are separated by a finite gap from the energies of the eigenstates that do not belong to M. The projector P is the 
projector onto the model space M. Following the notation of the Introduction, the energies of the eigenstates that 
belong (resp. do not belong) to M are denoted by Ef (resp. Ef). 

In this section, we prove the convergence of each term of the perturbation expansion of the wave operator when 
e — > 0. For notational convenience, we assume that t < 0. We first give a nonperturbative proof of this convergence. 
Then, we expand in series and we consider the different cases of the previous sections. 

A. Nonperturbative proof 

The nonperturbative proof is important in this context because its range of validity is wider than the series expansion 
(no convergence criterion for the series is required). Moreover, the proof that the wave operator indeed leads to an 
effective Hamiltonian is much easier to give in the nonperturbative setting. 

The first condition required in the nonperturbative setting is that the perturbation V must be relatively bounded 
with respect to Hq with a bound strictly smaller than 1. This condition is satisfied for the Hamiltonian describing 
nuclei and electrons interacting through a Coulomb potential 3 ^. Before stating the second condition, we define 
the time independent Hamiltonian h(X) = Hq + XV, its eigenvalues Ej(X) and its eigenprojectors Pj{X), such that 
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h(X)Pj(X) — Ej(X)Pj(X). The second condition is that the eigenvalues Ej(X) coming from the eigenstates of the 
model space (i.e. such that Pj(0)P = Pj(0)) are separated by a finite gap from the rest of the spectrum. According 
to Kato^, the eigenvalues and eigenprojectors can be chosen analytic in A. Then, a recent version of the adiabatic 
theore m 39 ' 40 shows that there exists a unitary operator A, independent of e, such that 

lim||£/ e (0,-oo)Pj(0) -e tB ^ e APj (0)| = 0, 
where U e (t,to) is the evolution operator for the Hamiltonian H e (t) — er t \ t \Q bH ° t Ver %Hat and 

In other words, the singularity of U e (0, — oo)Pj(O) is entirely described by the factor e l6 ^ e . The operator A satisfies 
the intertwining property AP/(0) = Pj(l)A. To connect this result with the case that we are considering in this 
paper, we have to choose a model space M that satisfies the following condition: there is a set I of indices j such that 
P = J2jei -Pj(0)j where P is the projector onto M. 

This enables us to give a more precise condition for the invertibility of PC/ e (0, — oo)P. The adiabatic theorem 
shows that, for small enough e, PU e (0, — oo)P is invertible if and only if PAP is invertible. If we rewrite PAP = 
J2j PAPj{Q) = J2j p Pj(^) A : tne unitarity of A implies that PAP is invertible iff the kernel of £\ PP 3 {1) is trivial. 
We recover the well-known invertibility condition 1 ^ that no state of the model space should be orthogonal to the vector 
space spanned by all the eigenstates of H with energy Ej(l), where j runs over /. Note that, when the condition of 
invertibility is not satisfied, it can be recovered by adding the perturbation step by step^ 9 -. 

Then, we have 

Theorem 6 With the given conditions, the wave operator 

H := limC/ e (0,-oo)P(PC/ e (0,-oo)P)- 1 

is well defined. Moreover, there are states \(pj) in the model space such that Q\(pj) is an eigenstate of H with eigenvalue 
Ej(l) and the effective Hamiltonian H e s := PHQ satisfies H c g\(pj) — Ej(\)\(pj) . 

Proof. We first define Aj k = Pj(0)AP k (0), for j and k in /. Then an inverse B of PAP is defined by J^kei Aj k B k i = 
SflPjiO) where B jk = ^(0)5^(0). Then, (PU e (0, -co)/ 3 )- 1 ~ £ Jfc e~^^B jk and 

fi £ (0,-oo) := [/ e (0,-oo)P(P[/ e (0,-oo)P)- 1 ~^APj{0)B jk . 

Since the right hand side does not depend on e, then O e (0, — oo) has no singularity at e = and 

Q. = limfi e (0,-oo) = ^APj(0)B jk . 

This proves the existence of the wave operator. To prove the existence of the states \(p 3 ) of the theorem, define 
\<fij) — PA\<fj), where \<pj) is an eigenstate of Pj(0): Pj(0)\(fj) — \<pj). Indeed, we have 

U\<Pj) = UPAPjiO)^) = ^£ / AP k (Q)B km A mj \<p j ) = AR,(0)|^) = ^(1)^1^), 

km 

where we used the intertwining property in the last equation. We can now check that f2|y>j) is an eigenstate of H 
with eigenvalue Ej(l). 

HQ\^} = = ^(1^(1)^) - ^(lp^). (9) 

Finally, by multiplying eq. Q by P from the left, we obtain 

ffeffl^) = EjMPUlw) =^(1)|^-), 

because Pfl — P and P\<fij) — \<Pj)- n 

Thus, the eigenvalues of H e Q are eigenvalues of the full Hamiltonian H . This is exactly what is expected from an 
effective Hamiltonian. In practice, the operator A is not known and the states \ipj) are obtained by diagonalizing 

//,::• 
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B. Series expansion 

We consider again the series expansion in terms of permutations. A straightforward calculatio n 41 ' 42 of the Picard- 
Dyson series gives us 

TT(n nuf r \^ 1 )(^ 1 \v\^)---{^- 1 \v\^J(^\v\<s>°) 

C/e(0,-Oc)|$ > = |$ 0> + Ijr _ + _ ^ + {n _ _ {Eo _ ^ - ~ ) , 

where we used the completeness relation 1 = This expression clearly shows that the terms of the 

expansion (and the evolution operator) are divergent as e — > when any Ei k is equal to Eq. 
For (j € S n , we set fl a (t) :— ft a (t, — oo). We then have 

n a (t) = {-i) n j dti / 1 dt 2 ... / " 1 dt„Qe (e+lffo)t "< 1 >V r e- lffot -( 1 'i? 1 

J —oo J —oo J — oo 

e (e+iH )t a{2) y e -~iH t a{2) jj2 _ ^ ^n-l e (£+iH )t <T („)^/ e -2-ffot <T („)p j 

where i?£ := Q if <r(fc + 1) > cr(fc) and J?J := -P if er(fc + 1) < a(k). We replace i?£ by ±^ Qfc+1 where, 

if R% = Q, then ± = + and the sum is over the image of Q, and if i?£ = — P, then ± = — and the sum is over the 
image of P. Thus 

n a (t) = {-i) n (-l) d j dt x j 1 dt 2 ... j " 1 dt n e {e+lFl ^ F2)t ^ 1 ->e^ +lF2 - lF ^ t ^ _ _ e (t+iF n -tF n+1 )t a{n) 

J —oo J — oo J — oo 

«2) ■ • ■ (a n \V\a n+1 )(a n+ i\, (10) 

ai...a„+i 

where d is the number of elements of the descent set of a, Fi is the energy of a,i and where the sum over a,\ is over 
the image of Q, the sum over a n +i is over the image of P and the sum over a k for 1 < k < n + 1 is over the image of 
Q if er(fc) > a{k — 1) and over the image of P otherwise. Consider now the time integral 

f a {t) := /' dt x F dt 2 ... /"'"^dt^+^-^^Me^^^ 



-oo J —oo 

t 



, (1) /" 4 T(n-l) 

dftr(i) / d Sr(2) . . . / d^e^^^e^ 2 -^ 2 . . . e^+^-^H 



— oo 



where t = a 1 . The integral over s T („) is 

7-oo " " ' ( e + ^Fr{n) - *F T („) + i) 

The integrand of the integral over s T („_i) becomes 

e (2e+i(F T( „ ) +F T(n _ 1) -F T(n)+1 -F T („_ 1)+1 )s T( „_ 1) 

(e + «F r(n) - iF T (n) +1 ) 
A straightforward recursive argument shows that 



x a {\) . . . x a { n y 

where 

X a (k) := ke + i(F a -i^ n) H h F„-\^ n _ k+1 ) ~ F<r-i(n)+i ^r-i(n-fc+i)+i)- (11) 

Therefore, 

«a(t) = E Vm l"i)("i|^l«2)...(a n |F|a n+1 )(a ra+1 l. (12) 
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C. Examples 

A few examples of tt a (0) are 

fi(i)(0) = ( _ t) !**><**! *W<*il 



(12)( } = ( "° *k (^?-m(^(E?-E D y 



n ( 2i)(o) = -H) 



Finally we consider two examples that will prove useful: 

o (Q) = _ ( _ t) 3 v m{<P4v\^){^\v\^k)^ k \v\^i)(^\ 

( " 3) ( e + < E % E H) (2e + *(^ Q + J5? - £f - Ef)) (3e + ,(J5? - 25f )) 



0(312) (0) - -( *) 3 ^ , ? jr/'.Wo. ■ r'<< ^ rra,,, . ,,,r< r 



. , ,, , , , i>( (e + - Ef )) (2e + ,(£« + ^ - £f - Ef)) (3e + - Eft) 
By adding these two terms, we obtain a denominator involving only the difference of two energies. 

\$i)(&i\v\$ j )($ j \v\$ k )(&k\v\$i)($i\ 

( e + < E ? ~ E D) ( e + < E k - Ef)) (36 + i{E? - Eft) ' 



Note that the sum is simpler than either ri( 2 i3)(0) or fi(3i2)(0). This is a general statement and the simplification 
becomes spectacular at higher orders. For the example of n = 7, there is a single tree T which is the sum of 80 
permutations, and the denominator of Qt is simpler than the denominator of VL a for any of the 80 permutations a of 
St- This will be proved in section ITVl Note also that, if we assume that the states in the image of P (i.e. the model 



space) are separated from the states in the image of Q by a finite gap 5, so that Ef — Ef > 6, then the denominators 
of all the examples are non-zero when e — > 0. In other words, the limit lim^o fi CT (0) exists for all the examples. In 
the next section, we show that this result is true for all permutations a. 



D. Convergence of £l a (t) 

Definition (|11[) is convenient for a computer implementation but it does not make it clear that X a (k) is nonzero 
if e = 0. For that purpose, we need an alternative expression for X a (k), which is essentially a corrected version of 
the graphical rule given by Michels and Suttorpi^. We first extend any permutation a £ S n to the sequence of n + 2 
integers a = (<x(l), • ■ • , o(n + 2)) = (0, <r(l), . . . , cr(n), 0). Then, for k S {1, ... , n}, we define the two sets 

S*<(fc) : = {i\l<i< n + l and a(i) < k < a(i + 1)}, 
S*>(fc) : = {i\l<i< n + l and a(i) > k > a{i + 1)}. 

For example, if a = (41325), then 5<(1) = {1},S<(2) = {1,3},5<(3) = {1, 3, 5}, 5<(4) = {1,5},5<(5) = {5} and 
S>(1) = {6},S>{2) = {2,6},S>(3) = {2,4,6},S>(4) = {2,6},£>(5) = {6}. The graphical meaning of these sets is 
illustrated in Figure [5] Notice that the vertical axis is oriented downwards in order to reflect the time-ordering in the 
integrals S a (t). 

Lemma 7 (i) S„(k) cannot be empty and (ii) S^r(k) and S^(k) have the same number of elements. 

The lemma follows from the graphical interpretation of the construction of S>(k) and S^(k). The graph of a 
(constructed as in Figure 1) is a sequence of edges connecting the points (i,a(i)). Since the graph starts from (0,0) 
and since there exists one point with ordinate n, any horizontal line with non integer ordinate y, < y < n, will be 
crossed from above by a segment (remember the vertical axis is oriented downwards) . A similar elementary topological 
argument shows that such a horizontal line is always crossed successively from above and below by segments, the 
series of crossings starting from above and ending from below, which implies |S<(fc)| = \S^(k)\. 

The key step in the proof of convergence is 
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1 2 3 4 5 6 7 


1 

2 
3 
4 
5 

FIG. 2: Construction of (3) for a — (41325). We build a — (0,4, 1,3,2,5,0), we draw a continuous line L starting from 
(1,ct(1)) = (1,0) to (2,ff(2)) = (2,4), to (3,a(3)) = (3,1), up to (7,cr(7)) = (7,0). To determine S<(3), we draw a 
horizontal dashed line just above 3 and we gather the segments of L crossing the dashed line from above. In our case the 
segments are ((1, 0), (2, 4)) , ((3, 1), (4, 3)) and ((5, 2), (6, 5)) . (3) is the set of abscissae of the first point of each segment: 
(3) = {1, 3, 5}. Similarly, (3) is obtained from the segments that cross the dashed line from below: (3) = {2, 4, 6}. 

Lemma 8 For a G S n and a given sequence a\, . . . , a„+i, compatible with a (see eq. klfljl ). we have 

X a (k) = ke+ ]T iF?- Y, lF f> ( 13 ) 

jeS<(n-k+l) jeS>(n-k+l) 

where we write F® — Fi (resp. Ff = Fi) when on belongs to the image of Q (resp. of P). 

Proof. We first show that it is true for k = 1. Indeed, X a (l) = e + i{F a -i^ — F a -ir n )+\)- Let j — cr _1 (n), we have 
o~(j) — n. Then, either j ' = 1 and \a\) is in the image of Q, or j > 1 and a G S n implies o~(j) = n > a{j — 1), so that 
\a.j) is in the image of Q. Thus, in all cases, F a - = F^-xi n y Consider now Either j = n and |aj+i) = |a n +i) 

is in the image of P, or j < n and <t(j + 1) < cr(j) = n, so that |o!j+i) is in the image of P. Thus, in all cases, 
F<r-i(n)+i = F a-i( n )+v Therefore, X„(l) = e + i(F®_ 1{n) - F p _ l(n)+1 ). On the other hand, S<(n) = {cr" 1 ^)} and 
S> (n) = {cr _1 (n) + 1} since (j, n) is the only point of the graph with ordinate n. Thus, the two members of eq. (|13j) 
arc equal for k = 1. 

Assume now that eq. (|13[) holds for all the X a (i) with i = 1, . . . , fc, k < n, and consider the equation (which is true 
by definition of the X (7 (i)s): 

X CT (fc + l) = X ff (fc) + e + *(F i -.F J - + i), (14) 

where j — cr _1 (n — k). We first treat the case 1 < j < n. Four possible situations can arise: (i) a(j — 1) > ct(j) > 
a(j + 1), (ii) a(j - 1) < <r(j) > a(j + 1), (hi) tr(j - 1) > cr(j) < <r(j + 1) and (iv) ct(j - 1) < cr(j) < a\j + 1). In 
case (i), we have Fj = Ff and Fj+i = Ff +1 . On the other hand, condition (i) implies a{j) > a(j + 1) > a(j + 2), so 
that S>(n — k) = S>(n — k + 1) and S^(n — k) is obtained from Sf{n — k + 1) by removing {j} and adding {j + 1}. 
eq. (fT4|) together with the hypothesis that eq. (fl"3"| holds for -Xo-(fc) imply that eq. (fT3|) holds for A CT (fc + 1). Case 
(ii) implies Fj = F9 and Fj+i = Ff +1 , case (iii) implies Fj = Ff and Fj+\ = Ff +1 , case (iv) implies Fj — Fj and 
Fj + i — Fj +l . In all cases, these identities imply that the two expressions (fT4"]) and (fT3")) for -Xo-(fc + 1) do agree. 

It remains to treat the boundary cases. If j = 1, then Fj = F® and we have either (i) <j(l) < <r(2) and F 2 = F® 
or (ii) a(l) > cr(2) and F 2 = Ff . We know that — 0, thus, case (i) corresponds to a(l) < ct(2) < <r(3), so that 
according to eq. (fl"3|) . 

A CT (fc + 1) - X CT (fc) = - F?), 

in agreement with eq. (fT4"|) . In case (ii) we have tr(l) < ct(2) > <r(3), which amounts to add iF® and remove 
Again, eq. (fl3|) holds for fc + 1. Finally, if j = n, then Fj+i = -P^+i and we have (i) a(n — 1) < a(n) and F n = F® 
or (ii) a(n — 1) > <r(n) and F„ = Ff. Case (i) corresponds to a(n) < a(n + 1) > (r(n + 2), case (ii) corresponds to 
a(n) > a(n + 1) > a(n + 2). In all cases, the relation given by eq. (fT3| is satisfied for k + 1 and the induction proof 
is complete. □ 

We can now ready to prove 
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Theorem 9 The limit 

Q a = limQ CT (0), (15) 

e-yO 

is well-defined. 

Proof. To prove the term-wise convergence of £7 CT (t) as e — ► 0, consider eq. (fl~3|) . The sets S<(n — k + 1) and 
S>(n — k + 1) have the same number of elements, say rik and, by the gap hypothesis, we have F- — -F 1 / 3 > (5 for any 
z and j. Therefore, |X CT (A;)| 2 > k 2 e 2 + ri^S 2 > S 2 , since > 1. Thus, the denominator remains away from zero by a 
finite amount for any e > and the limit of 1/X a (k) for e — > is well-defined. □ 



IV. TREES 



We showed that, for each permutation a, the wave operator fi CT (i, — oo) has a well-defined limit as e — > 0. The 
detailed proof was rather lengthy and the final expression for X a (k) suggests physically the simultaneous occurrence 
of transitions from states of the model space to states out of it. The convergence is actually much easier to show 
in terms of trees, and the expressions showing up in the expansion are simpler, mathematically and physically, each 
factor of the denominator corresponding to a single difference between an energy in the model space and an energy 
out of it. 

In this section, if N is the dimension of the model space, we write i € Q for i > N and j £ P for 1 < j < N, both 
for notational simplicity, and to emphasize the meaning of the indices, that correspond respectively to eigenstates in 
the image of Q and in the image of P (i.e. the model space). 

Proposition 10 IfT = Ti V T 2 , then, for t < 0, 

fi T (t):=fi T (t,-oo) = Y, e^-^+m^l^)^-!, (16) 
ieQ, 3 eP 

where Qj, is obtained recursively by: 
For T = Y , fi?. := -i .ltf^ 



Y 



ForT x = \ ,T 2 J= | ; ilf, := -l \, 



ForT 1 ^\,T 2 =\:n«:= l Z °g ^ ^ . 

I iE?-iEf + \T\e 



keP 



For Ti ^ | , T 2 ^ 



Proof. The computation of f2 i,? follows from eq. (fT2|) . Let us consider for example the case T\ = \ ,T%=£ [ . Then, 
by applying theorem [3] 

n T (t) = -i Y f dsQe es e- tE ? s e {lE ?- iE ? + \ T ^ s n%e* HoS V\$ t )($ 3 \. 
ieQjepJ-°° 

We replace Q by J2 l$fc)($fc| and obtain, by using \T\ = \T 2 \ + 1, 
fceQ 

n r (t) = -* X! / d se ^ Q -^ p +i T i e ) s ($ fc |i/|$ 4 )^J<i> fc )(^| 

k,ieQ,jeP^-°° 

(fBS-tBf + \T\e)t 

k,ieQ,jeP %bj k - tL j + I 1 \ e 

The two other cases can be treated similarly. □ 

Since, for arbitrary i and j, E 1 ^ — > 6, we get: 
Corollary 11 The limit Clx(t) = limfJiYO is well defined. 
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A. Relation with the Rayleigh-Schrodinger perturbation theory 

Kvasnicka 7 and Lindgren 6 independenty obtained an equation for the time-independent Rayleigh-Schrodinger per- 
turbation theory of possibly degenerate systems: 

[uJ,H a ]P = VujP - ZJPVuJP, (17) 

where the time-independent wave operator ZJP transforms eigenstates |$o) of Hq into eigenstates uJP|$o) of Hq + V 
and where PujP = P (see ref.|H p. 202 for details). 

Equation (jTT)) is an important generalization of Bloch's classical results^ because it is also valid for a quasi- 
degenerate model space (i.e. when the eigenstates of Hq in the model space have different energies). The relation 
between time-dependent and time-independent perturbation theory is established by the following proposition: 

Proposition 12 We have UP = f2 = lim e _j.o £1(0). 

Proof. We take the derivative of eq. (|5J) with respect to time and we substitute \(t,to) — Cl(t,to) — P- This gives us 

i^-fi(t,t ) = H(t)n(t,t )-si(t,t )H(t)n(t,to). 

at 

If we take t — and to = — oo, we obtain by continuity 

i-r^(t)\ t =o = VWO)-fi(O)VWO). 
at 

When we compare this equation with cq. (|17p . we sec that uJP and Q satisfy the same equation if 

zhm^| f=0 = [limn(0),flb]. (18) 

e^O at e— >0 

To show this, we prove it for each term £It- Indeed, eq. (TlB)) gives us 

j d^_(t) = ^ (Ef - + i \T\e)e^-^+\ T ^\^)(^l 
ieQ,jeP 

and 

[n T (t),H ] = Yl (Sf--E?)e^-< +|T|e ^>i)<*i|. 
ieQ,jeP 

By continuity in e, these two expressions are identical for all t when e —> 0. If we take t — and we sum over all trees 
T, then we recover eq. (fT5|). Therefore, uJP and satisfy the same equation. It remains to show that they have the 
same boundary conditions: fiP = and PCI = P. By eq. ([3]), these two equations are true for Q(t) with any value 
of t and e. □ 

As a corollary, proposition 1101 provides a recursive construction of the wave operator f2. 



B. An explicit formula for ilr 



In this section, we show how fir(i) can be obtained non-recursively from the knowledge of T. The key idea is to 
replace T by another combinatorial object, better suited to that particular computation. We write therefore 7t for 
the smallest permutation for the lexicographical ordering in St (we view a permutation as a word to make sense of the 
lexicographical ordering: to (35421) corresponds the word 35421, so that e.g. (35421) < (54231)). Since St is always 
non empty, the map 7 : T 1 — > jt is well-dehned and an injection from the set of trees to the set of permutations. 
These permutations are called Catalan permutations^, 312-avoiding permutations (i.e. permutations for which there 
does not exist i < j < k such that a(j) < <j(k) < cr(i), see ref. [3l| p. 224), Kempf elements^ or stack words^. 

Various elementary manipulations can be done to understand such a map. We list briefly some obvious properties 
and introduce some notation that will be useful in our forthcoming developments. If I — (ai, ...,cifc) is a sequence of 
integers, we write I[n] for the shifted sequence (a\ + n, a& + n). 
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Then, let T — T\ V T 2 be a tree. The permutation j(T) (that we identify with the corresponding word or sequence) 
can be constructed recursively as 7( | ) = 0, 7( Y ) = (1) and 

7 (T) := (7C7i)[l], l, 7 (r 3 )[|Tx| + 1])). 

The left inverse of 7, say T, is also easily described recursively as T(0) = | , T(l) := Y and 

7» := T((a(l), a(k))[—l)) V T((a(k + 2), ...,a(n))[-k 1]), 

where a € S n , a = (c(l), cr(fc), 1, a(k + 2), er(n)) is in the image of 7. 

Permutations in the image of 7 can be characterized recursively similarly: with the same notation as in the previous 
paragraph, a is in the image of 7 if and only if (c(l), a(k))[— 1] and (cr(fc + 2), a{n))[— k — 1] are in the image 
of 7, where k is the integer such that a(k) — 1, 

We are now in a position to compute f2-r(i)- Let us write At for all the sequences a = (oti, a n +i) associated to 
7(T) as in equation (|T0|) . We write, as usual, -Fj for the eigenvalue associated to cii. Recall that a± £ Q, a„+i G P 
whereas aj, i 7^ 1, n + 1 is in Q if <r(i) > er(« — 1) and in P otherwise. 

These sequences are actually common to the expansions of all the ^l a (t), a € St (this is because they depend only 
on the positions of descents in the permutations a £ St, as discussed in the proof of thm. [5]). They appear therefore 
in the expansion of Ot(0 = ^trif)- They actually also correspond exactly to the sequences of eigenvectors that 

ae S T 

show up in the recursive expansion of fi-r^) (proposition I16[) (this should be clear from our previous remarks, but 
can be checked directly from the definition of the recursive expansion) . We can refine the recursion of proposition 1101 
accordingly: 

Lemma 13 We have: flT(t) = ^e( lFl_lF "+ 1+ l T l e )*Q"|ai)(a n +i|, where fl"(t) is defined recursively by: 

a 

Y iE?-iE P + e 



4°"-' afc) (^|y|a fc+1 )47 +1 '-' a - a " +l) 

where T = T\ V T 2 and k = \Ti\. For T\ = \ , we have f4? = —1 and for T 2 — \ , we have Q^ n+1 ' = 1. 

Let us now consider jt- For i = 1, ...,n, we set: l(i) = inf{j < i\ik, j < k < i,7x(£;) > 7t(*)} and r(i) — sup{j > 
> I > i,7t(0 > 7t(*)}- In words, l(i) is defined as follows: consider all the consecutive positions k on the left 
of position i, such that the value of the permutation 7t(&) is larger than 7t(«)- Then, l(i) is the leftmost of these 
positions k. Similarly, r(i) is the rightmost position j such that, on all positions k between i and j, the permutation 
7t(&) is larger than 7t(«)- 

Theorem 14 We have: 



Q$ = (-t) n (-l) d - 1 <ai|V|a 2 )...(an|^|a n+ i)n- 



or, equivalently, 

*M*) = HH-l)^ 1 ^ \a 1 ){a 1 \V\a 2 ).4a n \V\a n+1 )(a n+1 \e^-^+^ t 

n 1 

niF l{:l) -iF r[l)+1 + (r(i)-l(i) + l)e 
where d is the number of leaves of T pointing to the right. 

For example, if T = J, then lT = (213), Z = (113), r = (133) and 

V*^*, ( e +^ Q -^))( e +^ fc Q -^ p ))(3 £ + i (^ Q -i?r))' 
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Proof. We show that fi^ as defined in theorem [14] satisfies the recursion relation of lemma [T3l Let us first consider 
that T\ ^ | and T% ^ | . Let iq denote the index such that 7t(*o) — L with 1 < io < n. The term of the 
numerator corresponding to io is (ai \V\a.i +i), which is the central term of the recursion relation. We have l(io) = 1 
and r(io) = n. Thus, the denominator is «-F/(i ) — iF r ^ +1 + (r(i ) — l{io) + l) e = tFi — iF n+ i + ne, which is the 

denominator of the recursion relation. Now we check that the product of terms for i < io in theorem[T?]is fi^ 1 ' "' Ql ° _1 \ 
The matrix elements (a± \V\0t2) ■■■{oti Q -\ \V\a.i a ) obviously agree, so we must check that the denominators agree. Thus, 
we verify that, for 1 < i < io, lr(i) = ^T t {i) and tt{i) = r^ii), where It and tt denote the I and r vectors for tree 
T. We know that 7t(«) = 7Ti(«) + 1 for 1 < i < io- Thus, for k < i, 7t(&) > 7t(«) if and only if 7Ti(fc) > 7Ti(*) and 
= ^Ti(*)- For rr, we notice that for I = io we have 7t(0 = 1 < 7t(«) and the relation 7t(0 > 7t(*) does not 
hold. Therefore, r<r(i) < io and rx(i) = ^Ti(*) by the same argument as = ^7\(*)- The same reasoning holds for 
T2 and the recursion relation is satisfied. The cases T\ = \ or T2 = \ are proved similarly. □ 

We conclude this section with a geometrical translation of the previous theorem. Consider a tree T with \T\ = n 
and number its leaves from 1 for the leftmost leave to n + 1 for the rightmost one. For each vertex v of T, take the 
subtree T v for which v is the root. In other words, T v is obtained by chopping the edge below v and considering the 
half-edge dangling from v as the dangling line of the root of T v . For each tree T v , build the pair (l v , r v ) which are the 
indices of the leftmost and rightmost leaves of T v . Recall that \a.i) belongs to the image of Q (resp. P) if leaf i points 
to the left (resp. right): this implies in particular that Fi v = F® and F Tv = F^ v - Then, the set of pairs (l v ,r v ) where 
v runs over the vertices of T is the same as the set of pairs (l(i),r(i) + 1) of the theorem. The formula for J7y can 
be rewritten accordingly and this geometrical version can be proved recursively as theorem 1141 Conversely, it can be 
used to determine the tree corresponding to a given denominator. 

V. CONCLUSION 

We considered three expansions of the wave operator and we proved their adiabatic convergence. We proposed to 
expand the wave operator over trees, and proved that this expansion reduced the number of terms of the expansion 
with respect to usual (tractable) ones, simplified the denominators of the expansion into a product of the difference 
of two energies and lead to powerful formulas and recursive computational methods. 

We then showed that this simplification is closely related to the algebraic structure of the linear span of permutations 
and of a certain convolution subalgebra of trees. 

As far as the many-body problem is concerned, we showed that the simplification of diagrams is not due to the 
details of the Hamiltonian but to the general structure of the wave operator. When the eigenstates and Hamiltonian 
are expressed in terms of creation and annihilation operators and quantum fields, the algebra of trees mixes with the 
Hopf algebraic structure of fields^^. It would be interesting to investigate the interplay of these algebraic structures. 

The terms of the Rayleigh-Schrodinger series are usually considered to be "quite complicated" (see ref.H^, p. 8) and 
difficult to work with. The general term of the Rayleigh-Schrodinger series for quasi-degenerate systems is obtained 
as the limit for e — > of ^t(O) in theorem 1141 Through our recursive and non-recursive expressions for these terms, 
many proofs of their properties become almost trivial. The tree structure suggests various resummations of this series, 
that will be explored in a forthcoming publication. 
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Appendix A: A crucial lemma 

In this appendix we state and prove a lemma that is crucial to demonstrate theorem [5l which is one of the main 
results of our paper. We first need a noncommutative analogue of Chen's formulas for products of iterated integrals. 
This analogue, proved in refs. [Ho| and [22|, provides a systematic link between the theory of iterated integrals, the 
combinatorics of descents, free Lie algebras and noncommutative symmetric functions (see ref. [22| and Appendix [Bl of 
the present article for further details). 
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Let L = (Li, L n ) be an arbitrary sequence of time-dependent operators Lj(i), satisfying the same regularity 
conditions as H(t) in section [TT1 Let a be a permutation in S n and define: 

n£(t,i ) := f dh f'dh... r 'dtnii (*„(!))...£„(*„(„)). 
•'to "to "to 

The notation is extended linearly to combinations of permutations, so that for /i := Y f 1 a-o' 1 with \x„ G C. 

Then H,^(t,to) is defined as the linear combination l^<r^(t,to). For K := (K±, K m ) another sequence of 

cr65„ 

time-dependent operators, we write L ■ K for the concatenation product (Li, L n , K\, K m ). 

We also need to define the convolution product of two permutations. If a G S n and (3 G S m , then a * /3 is the sum 
of the ("^ m ) permutations 7 G <S n +m such that st(j(l), 7(77.)) = (a(l), a(n)) and st(^(n + 1), 7(71 + m)) = 
(/3(1), ...,/?(m)). Here, st is the standardization map defined in section Hi C II For instance, 

(2,3,1)*(1) = (2, 3, 1,4) + (2, 4, 1,3) + (3, 4, 1,2) + (3, 4,2,1), 

(1,2)* (2,1) = (1,2,4,3) + (1,3,4,2) + (1,4,3,2) + (2,3,4,1) + (2,4,3,1) + (3,4,2,1). 

In words, the product of two permutations a G S n and j3 G S m is the sum of all permutations a of S n+m such that the 
elements of the sequence (cr(l), . . . , cr(n)) are ordered as the elements of (a(l), . . . , a(n)), in the sense that a(i) > a(j) 
if and only if a(i) > a(j) and the elements of (<r(n+ 1), . . . , a(n + m)) are ordered as the elements of (/3(1), ■ ■ • , /3(m)). 
Now, we can state the noncommutative Chen formula^ (see also remark 3.3, p. 4111 of ref. 1221 ) 

Lemma 15 We have: 

Q£(Mo)Of(t,t )=fl^(Mo)- 
The following lemma can be proven similarly: 

Lemma 16 We have, for L and K as above and J a time- dependent operator: 



dstl^s, t )J(s)ilf(s, t ) = J2 n£- (J) -*(t, to), 

7 



where 7 runs over the permutations in iS n + m +i with 7(71 + 1) = 1, si(7(l), ^(n)) = a, st(^f(n + 2), 7(71 + 771 + 1)) 
3- 

Proof. We expand fi^ and Vt^ 



dsCt a (s,to)J(s)Qp (s, t a ) = / ds / dwi . . . / dw„ / dvi . . . dv m 

to «to •'to «to "to •'to 

■ • ■Ln(w Q (n))^(s)ii'l(^(l)) ■ ■ .if m (t>8( m )). 

This is the same formula as for the expansion of ^^(s, to)D,-p(s, to), up to the term J(s) and the integration f t ds 
that however do not change the underlying combinatorics. Therefore, lemma [15] holds in the form 

ft pt rs n + m -i 

dsn^(s,t )J(s)nf(s,to) = J2 ds dsi.../ ds n+m 

'to cr J to Jto 

n+m) J ; 

where, by the definition of a * /3, the sum over <r is over all the permutations of iS rl + rn such that st(a(l), . . . , o~(n)) = a 
and st(a(n + 1), . . . ,o~(n + to)) = f3. Now, we change variables to t\ = s, = Si for i = 1, . . . ,n + m. The 
permutation 7 of ti, . . . , /j n+m+1 corresponding to the permutation a of s%, . . . , s n + m is characterized by 7(72 + 1) = 1 
(because s > Si for all i), 7(1) = cr(i) + 1 for 1 < i < n and j(i + 1) = a(i) + 1 for ri + 1 < i < n + to. Therefore, 



/' 

J to 



t /-t /-tn 



dsn^(s,to)J{s)nf(s,t ) = 2^ / dti... / dt 



tO ^t0 ^t0 



^l(^7(l)) • • • ^ri(*7(n))^(*7(n+l))-^l(*7(n+2)) • • • ^m(t 7 (n+m+l))i 

where 7 satisfies 7(71 + 1) = 1, st(7(l), 7(n)) = a, st(^(n + 2), 7(77, + to + 1)) = (3. The lemma is proved. □ 

Notice that the same process would allow to derive combinatorial formulas for arbitrary products of iterated integrals 
and for integrals with integrands involving iterated integrals, provided these products and integrands have expressions 
similar to the ones considered in the two lemmas. 
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Appendix B: The algebraic structure of tree-shaped iterated integrals 

Our exposition of tree-parametrized time-dependent perturbation theory has focussed on the derivation of conver- 
gence results and explicit formulas for the time-dependent wave operator. However, the reasons why such an approach 
is possible and efficient are grounded into various algebraic and combinatorial properties of trees, descents and similar 
objects. 

These properties suggest that the theory of effective Hamiltonians is grounded into a new "Lie theory" generalizing 
the usual Lie theory (or, more precisely, generalizing the part of the classical Lie theory that is relevant to the study 
of the solutions of differential equations) . First indications that such a theory exists were already pointed out in our 
ref. HH. Indeed, we showed in this article that descent algebras of hyperoctahedral groups and generalizations thereof 
are relevant to the time-dependent perturbation theory. Applications included an extension of the Magnus expansion 
for the time-dependent wave operator. Our results below provide complementary insights on the subject and further 
evidence that algebraic structures underly many-body theories. 



1. The algebra structure 

We know from ref. and Sect III Bl that the family of integrals Q^(t,to) is closed under the product. This 
closure property is reflected into the convolution product of permutations. This result is a natural noncommutative 
generalization of Chen's formula for the product of iterated integrals. We show, in the present section, that the same 
result holds for integrals parametrized by trees. We will explain, in the next sections, why such a result -which may 
seem surprising from the analytical point of view- could be expected from the modern theory of combinatorial Hopf 
algebras. 

For L = (L%, L n ) a family of time-dependent operators (with the usual regularity conditions), and T a tree with 
n internal vertices, we write 

/* rH rt n -x 

dh dt 2 ... dt n L 1 {t a{1) )L 2 (t ry{2) )...L n {t a{n) ). 

,7 -OT ~ J 

This notation is extended to linear combinations of trees, so that e.g. if Z = T + 2T', where T and T' are two arbitrary 
trees with the same number of vertices, then f2§(t, to) = fi^(£, to) + 2f2j.,(i, to). For i < n, we write L<i = (L%, Zj), 

L>i = {Li, Ln). 

Proposition 17 For L = (L\, L n ) and K = (Ki, K m ) two families of time- dependent operators and T — T\\/T 2 , 
U = Ui V U2 two trees, \T\ = n, \T\\ = p, \T 2 \ = q, \U\ = m, \U\\ = I, \U 2 \ — k, we have: 

t t 

s%(t,to)n$(t,to) = J d s n{f{s,to)n^ l { s ,to)K l+1 ( s )n^ l+2 (s,t ) + J d s n^ p (s,t Q )L p+1 ( s )n^ p+2 (s,t )n«( s ,t ). 

to to 

In the formula, one or several of the trees Ti 1 T 2l U\,U 2 may be the trivial tree \ . 

Proof. Recall the integration by parts formula. For any integrable functions / and g we define F(t) := dsf(s), 
G(t) := dsg(s) and H(t) := F(t)G(t). Then 

H(t) = [ ds^^-= f dsf(s)G(s)+ f dsF(s)g(s). 

J to as Jto Jt 

Now, we use this identity with F(t) = f2j.(f , to). It follows from the proof of theorem [5] that f(s) is given by 

f(s) = Sl^(s,t )L p+1 (s)n^ +2 (s,to), 

with a similar identity for G(t) = fly (t, to). The proposition follows. □ 

In particular, it is a consequence of the proposition and a straightforward recursion argument that the linear span 
of the integrals Qy(t, to) is closed under products. 
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2. Hopf algebras and Lie theory 

This result may be formalized algebraically. Let us write T for the set of formal power series with complex 
coefficients over the set of trees. Proposition [T71 enables us to define a product on trees, denoted by *, such that 
f2j.(f, t )n^ (t, to) = Q,T*u(t'to)- This product is defined recursively by the equation 3 ^ 

T * U := (T * Ui) V U 2 + T x V (T 2 * U). 

Proof. The empty tree | is the unit for the product *. Assume that T* U is defined and satisfies the recursive relation 
for all trees such that \T\ + \U\ < n, and consider two trees T and U with \T\ + \U\ = n. The first term on the right 

hand side of proposition [IT] is J t * daft T (a,to)ft u ^ 1 (s,to)Ki + i(s)fl^ l+2 (s,to). By the recursive relation, we have 

n T (s,t )tt^'(s,t Q ) = il T '*^f. 

Thus, the whole term can be written 

Jt 

The second term of proposition [T7J is treated similarly and we obtain 

fl T (t,to)Fl§(t,to) — Q(t*Ui)VU2 + ^(T%{T 2 *U) 

Therefore, the relation Q T (t, t )^l^ (t, t ) = 0^y(t, i ) gives us 

T * U = (T * Ui) V U 2 + Tx V (T 2 * U). 

□ 

Corollary 18 The product provides T with the structure of an associative algebra. 

The corollary is a by-product of the associativity of the product of operators and of proposition 1171 

Of course, although the analysis of tree-shaped iterated integrals leads to a straightforward proof, the associativity 
propert y is a purely combinatorial phenomenon that originates ultimately from the associativity of the shuffle product. 
See ref. [22j for the connections between the noncommutative Chen formula and shuffle products, see also section 4 of 
ref. [m] for Schiitzenberger's classical (but rarely quoted) analysis of the formal properties of the shuffle product -in 
fact, the splitting of the convolution product of trees reflects the classical splitting of shuffle products into left and 
right half-shuffle products that had appeared in the study of Lie polynomials and was first encoded combinatorially in 
ref. The associativity can also be checked directly or deduced from the associativity of the convolution product * 
of permutations, since one may verify easily, using e.g. our description of the permutations in St that the convolution 
product as introduced above is nothing but the restriction to T of the convolution product on the algebra S = 

n c[s n }. 

This fact that the linear span of trees (often written PBT) defines a subalgebra of S for the convolution product 
(and even a Hopf subalgebra, referred to as the Hopf algebra of planar binary trees, whereas S is referred to as the 
Malvenuto-Reutenauer or Hopf algebra of free quasi-symmetric functions in the litterature) is well-known. It was first 
observed in ref. and further investigated in a series of paper a 34 i 36 ' 37 i 52 . There is however a slight subtelty here, 
since the embedding of PBT in S considered e.g. in ref. [30 is not the one we consider (another projection map from 
permutations to trees is used) so that PBT and T, although isomorphic as algebras (the product rule in PBT is the 
same as the one in T, see e.g. proposition 3.2 of ref. [3(1 where PBT is written CfY^]) do not agree as subalgebras 
of S. 

However, a corollary of this isomorphism is that the structures existing on PBT carry over to the analysis of the 
algebraic properties of tree-shaped iterated integrals. The existence of a Hopf algebra structure seems particularly 
meaningful since the classical applications of the theory of free Lie algebras to the analysis of differential equations 
can be rewritten using the framework of Hopf algebras (see e.g. the accounts in refs. [53II53 ). 

The same observation holds for the direct sum of the hyperoctahedral group algebras, that also carries naturally 
a Hopf algebra structur o 55 ' 56 , and which applications to time-dependent perturbative Hamiltonians were studied in 
ref. [22|. We leave for further research the investigation of the possible role of these Hopf algebra structures in the 
context of many-body theories. 
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3. Permutations, trees and descents 



Another meaningful observation, along the same lines, is that the three expansions we derived, based respectively on 
sequences of P and Q (as in the first expression for f2(t, t a )), on trees (third expression) and on permutations (second 
expression), reflect at the analytical level the existence of projection maps and embeddings between hypercubes (or 
the descent algebra), planar binary trees, and permutations. Following Viennolj^, these maps have been at the origin 
of modern enumerative combinatorics. We refer to refs. 30.34.36.37 for a detailed study of these maps that emphasizes 
the existence of underlying geometrical structures that go beyond the Hopf algebraic ones. 
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